The two-impurity Kondo model with spin-orbit interactions 
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We study the two-impurity Kondo model (TIKM) in two dimensions with spin-orbit coupled 
conduction electrons. In the first part of the paper we analyze how spin-orbit interactions of Rashba 
as well as Dresselhaus type influence the Kondo and RKKY interactions in the TIKM, generalizing 
results obtained by H. Imamura et al. (2004) and J. Malecki (2007). Using our findings we then 
explore the effect from spin-orbit interactions on the non-Fermi liquid quantum critical transition 
between the RKKY-singlet and Kondo-screened RKKY-triplet states. We argue that spin-orbit 
interactions under certain conditions produce a line of critical points exhibiting the same leading 
scaling behavior as that of the ordinary TIKM. In the second part of the paper we shift focus and turn 
to the question of how spin-orbit interactions affect the entanglement between two localized RKKY- 
coupled spins in the parameter regime where the competition from the direct Kondo interaction 
can be neglected. Using data for a device with two spinful quantum dots patterned in a gated 
InAs heterostructure we show that a gate-controlled spin-orbit interaction may drive a maximally 
entangled state to one with vanishing entanglement, or vice versa (as measured by the concurrence). 
This has important implications for proposals using RKKY interactions for nonlocal control of qubit 
entanglement in semiconductor heterostructures. 

PACS numbers: 71.10.Hf, 73.21. La, 75.30.Hx, 03.67.Mn 



I. INTRODUCTION 

The study of few-electron quantum dots, or so-called 
artificial atoms 1 , has grown into a major field of research, 
fueled by the promise of future technological applica- 
tions, as well as by problems in fundamental physics. 
Unlike an ordinary atom, the properties of a quantum 
dot can be manipulated by electrical gates, allowing for 
well-controlled studies of a whole spectrum of single- 
and many-particle phenomena, including one— and two- 
channel 3 Kondo effects, Coulomb blockade^, and Pauli 
spin blockade effects^. Recently, nanoscale devices with 
two spinful quantum dots connected to a conducting re- 
gion have become the focus of interest, following the 
breakthrough experiment by Craig et al£ (see also Refs. 

When the dots carry spin-1/2 magnetic moments, 
the problem is that of the two-impurity Kondo model 
(TIKM)2. 

In this model, two localized spins S12 of magnitude 
5 = 1/2 are coupled to a sea of conduction electrons via 
a spin exchange, 



H — -ffkin + J Si ■ <T 1 + J S2 ■ (T2- 



(D 



Here <Ji = tp^ (x^T^tp p{x{) are the electronic spin den- 
sities at the sites of the localized spins (with t the 
vector of Pauli matrices), J is the spin exchange cou- 
pling, and -ffkin is the kinetic energy of the electrons. 
To second order in J, the conduction electrons mediate 
a spin exchange interaction between the two localized 
spins, the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 



interactio: 
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-Hrkky — K(R)S\ ■ S2 



(2) 
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where the strength and the sign of the coupling K(R) de- 
pends on the distance R between the spins. The competi- 
tion between the RKKY interaction in ^ and the direct 
Kondo spin exchange in ([1]) is governed by the ratio of 
K(R) to the Kondo temperature Tr- ~ £>exp(— l/2ppj), 
where pf is the single-electron density of states at the 
Fermi level. When \K(R)\^> T K , the RKKY interac- 
tion dominates the direct Kondo exchange and will lock 
the two impurity spins into a singlet (for K(R) > 0) 
or a triplet state (for K(R) < 0). In the experiment by 
Craig et al£, the system consisted of two spinful quantum 
dots in a gated GaAs/AlGaAs heterostructure, coupled 
to an open conducting region. When tuning the coupling 
K(R) between the dots, the Kondo signature of the con- 
ductance across one of the dots was seen to vanish, or 
become strongly suppressed, indicating that the RKKY 
interaction between the dots dominates the Kondo inter- 
action, forming either a singlet state with no Kondo ef- 
fect or a screened triplet state with a much weaker Kondo 
effect^. 

When the strength of the RKKY interaction be- 
comes comparable to the Kondo temperature, \K(R)\ ~ 
Tk, there is a crossover between the antifcrromagnetic 
RKKY regime ("local singlet") and the regime where 
the triplet impurity state is Kondo screened ("Kondo 
singlet") 13 . If the system possesses particle- hole symme- 
try, the crossover is expected to sharpen into a second- 
order phase transition, controlled by a non-Fermi-liquid 
fixed point^. However, since this symmetry requires a 
high degree of fine tuning, the possibility to observe the 
corresponding quantum critical state was for a long time 
judged as rather unrealistic. However, in recent work by 
Zarand et alr^ it was shown that the critical state can 
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be stabilized against particle-hole symmetry breaking by 
using a device where the two dots are connected to two 
separate leads and RKKY-coupled via a magnetic insu- 
lator. 

In a variation of the proposal by Zarand et al^, we 
here consider a device scheme where the RKKY-coupling 
is mediated by a Coulomb blockaded auxiliary electron 
reservoir (see FIG. [2]). This is the same type of device 
that was studied in Ref. Ql| with focus on effects from 
charge fluctuations between leads and dots. This alter- 
native type of setup also allows for the study of effects 
from spin- orbit interactions on the critical behavior of the 
TIKM. This is the main theme of the present paper. As 
we shall see, the presence of spin-orbit interactions adds 
a new twist to the problem, with some quite interesting 
repercussions. Specifically, we shall argue that spin-orbit 
effects under certain conditions produce a line of critical 
points, exhibiting the same scaling behavior as that of 
the critical TIKM without spin-orbit interactions (up to 
RG irrelevant scaling corrections). 

The role of spin-orbit interactions in the TIKM is also 
of interest away from criticality, considering proposals for 
exploiting devices with RKKY-coupled quantum dots in 
spintronics, with possible applications for future quan- 
tum computing. In its most basic implementation, a 
quantum dot carrying spin-1/2 is here used to repre- 
sent a qubit, built from its two spin states ||), || ) 17 ' 18 . 
A key problem is that of connecting the qubits in such 
a way that their coupling and entanglement become 
controllable 19 . In the original proposal by Loss and 
DiVincenzo 1 ^, the spin qubits are coupled by a Heisen- 
berg exchange, produced by a tunable electrostatic bar- 
rier between neighboring dots. For the purpose of com- 
putation, however, one must be able to achieve fast and 
efficient control also of the coupling between electron 
spins on distant quantum dots. A major step towards 
this goal was taken by Craig et al& in their experimen- 
tal realization of the TIKM. As shown in the experi- 
ment, the coupling K(R) can be turned on and off by 
the electrical gates that control the energy levels of the 
dots, thus making possible a realization of a nonlocal 
two-qubit logic gate 2 * 3 . While there are other competing 
schemes for achieving nonlocal coupling of spin qubits, 
based on optical 21 or magnetic^ 2 - control, the RKKY- 
mediated two-qubit gate has its distinctive advantage in 
being a simple and easily scalable implementation. An es- 
timate by Rikitake and Imamura 23 , shows the two-qubit 
decoherence time (set by the electronic environment that 
mediates the RKKY interaction) to be well within the 
bounds for efficient gate operations. 

In this context, a central question is to understand 
how robust the RKKY coupling is against competing in- 
teractions. In recent work by Cho and McKenzie 2 ^ it 
was shown that the entanglement between two RKKY- 
coupled spin qubits (as measured by the concurrence^-) 
needs a minimum non-zero antiferromagnetic correlation 
determined by the competition between the RKKY inter- 
action and the Kondo effect. In the second part of this 



paper we take the Cho-McKenzie analysis 24 (sec also Ref. 
|26j) a step further by including the additional effect from 
spin-orbit interactions. A spin-orbit interaction mixes 
spin and charge and is known to be an insidious source 
of decoherence in spin-qubit devices^.. At the same time, 
it is possible that one could in fact exploit spin-orbit in- 
teractions for coherent control of qubit interactions, as 
proposed recently by several groups 2 ^. We here address 
the separate issue of how the presence of spin-orbit in- 
teractions influence the RKKY coupling and the entan- 
glement between two spin qubits. As we shall see, while 
the effect may be dramatic, it can be compensated for 
by properly tuning the electrical gates that define the 
device, and does not per se obstruct the operation of an 
RKKY-mediated two-qubit gate. 



The rest of the paper is organized as follows: In Sec. II 
we review how to derive a modified RKKY interaction in 
the presence of a spin-orbit interaction of Rashba type 29 . 
The analysis can easily be extended to a spin-orbit inter- 
action of Dresselhaus type, and we carry it out in parallel. 
The magnitude and sign of the modified RKKY interac- 
tion is seen to be strongly sensitive to the strength of the 
Rashba and Dresselhaus interactions, as encoded in the 
so called a- and /3-cocfficicnts 30 , with important implica- 
tions for the design of a spin-qubit device. In Sec. Ill we 
then turn to a study of how the Kondo interaction is influ- 
enced by spin-orbit interactions of both Rashba and Dres- 
selhaus type, capitalizing on recent work by Malecki 31 . 
In Sec. IV we use our results from the previous sec- 
tions to study the effect of the spin-orbit interaction on 
the quantum phase transition between the RKKY ("lo- 
cal singlet") and Kondo screened regimes. In the case 
when charge fluctuations are present we show that the 
spin-orbit interaction generates a line of fixed points (in 
the language of the renormalization group), and we pro- 
vide arguments for its interpretation. We also argue — 
on basis of symmetry arguments — that spin-orbit effects 
under certain conditions produce a line of critical points 
(unrelated to the line of fixed points just mentioned), pro- 
ducing the same leading scaling behavior as that of the 
critical TIKM without spin-orbit interactions. In Sec. V 
we change focus and explore how spin-orbit interactions 
influence the entanglement of two RKKY-coupled spin 
qubits (in the parameter regime where the direct Kondo 
exchange can be neglected). As a preamble we derive a 
general expression for the reduced density matrix and the 
concurrence for two spin qubits with U(l) symmetry, as 
appropriate for the present problem where the addition of 
spin-orbit interactions breaks the SU(2) symmetry of the 
standard TIKM down to U(l). Sec. VI contains a brief 
analysis of the behavior of the two-qubit concurrence at 
criticality, and in Sec. VII, finally, we summarize our 
results. 
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II. RKKY INTERACTION IN THE PRESENCE 
OF SPIN-ORBIT INTERACTIONS 



In the following we consider a realization of the TIKM 
in a nanoscale device where two spinful quantum dots are 
connected via tunnel junctions to a large central electron 
reservoir. The proposed device, which may be manufac- 
tured using a gated semiconducting heterostructure, is 
depicted in FIG. [I] The two dots are coupled to the reser- 
voir via point contacts, allowing the electrons to tunnel 
between dot i and the reservoir with amplitude Va,%- The 
dots are operated in the Coulomb blockade regime where 
transfer of charge between the dots and the reservoir is 
strongly suppressed but virtual fluctuations give rise to 
a spin-exchange (Kondo) interaction between the elec- 
trons trapped on the dots and the conduction electrons 
in the reservoir, as described by the Hamiltonian |T]). 
With Va,i = Va,2, the Kondo couplings to the dots are 
equal and given by J ~ V\ 1 /U, where U is the Coulomb 
blockade energy of the reservoir. As discussed in the 
introduction, there exists a parameter regime where the 
antifcrromagnetic RKKY interaction in © , generated by 
second-order Kondo exchanges, becomes dominant. As 
shown in Ref. l24l . a strong enough antiferromagnetic in- 
teraction leads to a finite (in fact, maximal) entangle- 
ment between qubits (spins on the dots) as required for 
a working quantum information device. 



interaction 



Vk- 
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i 

Dot 1 



Central 
reservoir 



Dot 2 



FIG. 1: (color online) The physical system that we study 
in this section. The different V are tunneling rates. The 
dots are operated in the Coulomb blockade regime, where 
charge transfer between the dots and the central reservoir are 
suppressed. 



We here focus on effects on the RKKY-interaction from 
the spin-orbit coupled electrons in the large central reser- 
voir. One may also enquire about the effect from the spin- 
orbit coupling of the electrons that reside in the quantum 
dots 1 and 2 (cf. FIG.[T]). While an interesting issue as 
such, we shall here bypass it by considering the idealized 
limit of two ultra-small spinful quantum dots, with the 
trapped electrons modeled as two completely localized 
spin- 1/2 impurities. In this limit only the spin-orbit cou- 
pling of the itinerant electrons in the large reservoir is to 
be taken into account. 

Spin-orbit interactions in semiconductor heterostruc- 
tures come in two guises, the Dresselhaus 32 and Rashba 33 
interactions, originating from the inversion asymmetry of 
the potential V(r) = V cr (r) + V'(r), where V cr (r) is the 
periodic crystal potential, and V'(r) is the potential due 
to confinement, impurities, external electric fields, etc. 
The electric field — VV(r) produces a Pauli spin-orbit 



H SO ~(px W(r)) • er 



(3) 



that can be of sizable magnitude due to the large po- 
tential gradients of the atomic cores 30 . Here er = 
(a x , <t v , <j z ) is the vector of Pauli matrices. If V(r) lacks 
inversion symmetry, i.e. V(— r) ^ V(r), then the Pauli 
interaction in f3| fails to average to zero in a unit cell, 
and results in a spin-splitting of the electronic bands. 
In semiconductors with bulk-inversion asymmetry (zinc- 
blende structures, including GaAs) the spin splitting can 
be encoded in the effective Dresselhaus interaction 32 



B k ^ k 



k})a\ 



(4) 



where ki,kj, and kg are the electronic wave numbers 
along the principal crystal axes, with cyclic per- 

mutations of (x,y, z), and where B is a material de- 
pendent coupling constant. For a heterostructure grown 
along [001], with the electrons confined to the a;y-plane, 
the Dresselhaus interaction reduces to 



H = (3(k x a x 



k y a y 



(5) 



when taking the average of Eq. ((4| along the z-direction: 
(k z ) — 0, and (fc 2 ) ~ (7r/d) 2 with d the characteristic 
electron wave length in the z-direction. Here f3 is the 
Dresselhaus /3-coefficient, j3 = —B(n/d) 2 . 

In a heterostructure, the spin degeneracy can be lifted 
also because of a structure inversion asymmetry of the 
confining potential contained in V'(r). This potential 
may obtain contributions also from an externally applied 
potential as well as from the effective potential from the 
position-dependent band edges. Assuming that — VV'(r) 
is an electric field along the z-direction, one obtains from 
([3]) the Rashba interaction 3 - 3 . 



a(k x a v - k y a x ), 



(6) 



where the Rashba coefficient a can be tuned via an ex- 
ternal gate^MS,. It may be worth pointing out that the 
spin-orbit interaction comes in two distinct varieties also 
in other materials. For example, in graphene, a Rashba 
interaction © controllable via an external gate electric 
field coexists with an intrinsic spin-orbit interaction de- 
termined by the symmetry properties of the honeycomb 
lattice, similar to the 2D Dresselhaus interaction ([5]) in a 
semiconductor heterostructure 3 -^. In the following we fo- 
cus entirely on effects produced by Hp and H a , as defined 
in Eqs. ([5]) and |6]) respectively. 

Turning to the RKKY interaction, it here acts between 
two localized spins in a two-dimensional (2D) elec- 
tron gas, where, in the device in FIG. 1, Si is attached to 
dot i (i = 1, 2), and the electron gas is represented by the 
central reservoir. It can be calculated in second order in 
the Kondo couplings J ii2 as^ 9 - 



H 



J1J2 



RKKY — 



Im / duj Tr [(Si ■ <t)G(R,uj + i0+) 

J — OC 



' —OO 

X 



S 2 -<x)G(-#,u; + iO + )], (7) 
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with u>f the Fermi level. Here G(R,u>) is the Green's 
function of a conduction electron with energy to, with 2R 
the distance between the localized spins (~ the width of 
the central reservoir in FIG. ([T])). The trace is taken over 
the two spin states of a conduction electron. The expres- 
sion in Eq. (0 can easily be adapted to the case when 
spin-orbit interactions are present by properly modify- 
ing the single-electron Greens function G. The case of 
a pure Rashba spin-orbit interaction was treated in Ref. 
129 . Here we will show that using a similar procedure as 
in Ref. 29, the form of the RKKY interaction can be cal- 
culated in the presence of both Rashba and Dresselhaus 
interactions of arbitrary relative strength. 

For this purpose it is convenient to write the 2D single- 
electron Hamiltonian with both kinds of spin-orbit inter- 
actions (cf. Eqs. ([H]) and ©) as 



H = 



2m 



-a\ k x 
a —pj \k v 



2m 



+ (Afe)-r, (8) 



where a and are the coupling strengths for Rashba and 
Dresselhaus spin-orbit interactions, respectivel y 32 ' 33 . All 
vectors that appear in Eq. 1(5)) have two components, i.e. 
fc = (k x ,k y ), t — (t x ,t v ). The scalar product of any 
vector m = (m x ,m y ) with the vector of Pauli matrices 
t is taken in the usual way, m ■ t = m x r x + m y T y . The 
Green's function corresponding to the Hamiltonian in Eq. 
(El) is 



G(k,u) = (w-iJ(fc)) -1 = {u 



2m 



(Ak)-T)- 1 



where 



= G (fc, lu; A) + Gi(k, oj; A) (Afc) • r, 

G (fc,w;A) = - ~~ 
Gi(k, to; A) = - 



(9) 



fc 2 > 

1m) 



(Afc) 5 



(10) 



h 2 



(Afc) 5 



Note that both Go and G\ are invariant under the trans- 
formation fc — * — fc. The Green's function in real space is 
obtained via a Fourier transform, where the spin depen- 
dence can be pulled in front of the integral by writing 



(Ak)e ik - R = -i(AV)e ltfi = I AR 
with R = R/\R\. We thus obtain 



d\R\ 



(11) 



G(R, to; A) = G (-R,w;A) + Gi(fi,w;A) (AR) -t, (12) 
where 

Gi(R, to; A) = -— / ,/-/,■ 



4tt 2 



e lk - R Gi(k,to;, 



d\R\ 



(13) 

Both Gq(R) and Gi(R) are invariant under the parity 
transformation R — > R. It is a straightforward task to 



perform the traces over the Pauli matrices in Eq. (0 (for 
details, see Appendix [A"]) . and by specifying a coordinate 
system where R = £ one obtains from Eqs. (0 - (|13jl the 
RKKY interaction in presence of Rashba and Dresselhaus 
spin-orbit interactions: 



rrSO _ tt 
"RKKY — -"Hois 



Hi 



Rashba 



H 



intcrf . : 



(14) 



where 



Hucis. 
-"Rashba 
"Dress, 
-^intcrf. 



F Si ■ S 2 

aFxiSx x S 2 ) v + a 2 F 2 S v S% 
0F 1 (S 1 x S 2 ) x +0 2 F 2 SfS 2 c 
a(3F 2 {SfS v 2 + SfS$) . 



(15) 



Here Ft = Fi(a, 0, R) are functions which in the gen- 
eral case are given by rather complicated integrals. For 
the case where only one type of spin-orbit interaction is 
present, they have been obtained analytically in the limit 
of large distances UfR 3> 1 and weak spin-orbit interac- 



tion a,0< §-fc 



One finds that 



F (0,a,R) =F (a,0,R) = 
Fi(0,a,R) =Fi(a,0,R) = 
F 2 (0,a,R) =F 2 (a,0,R) = 



F (0,a,R) 



sin2R\/k 2 F 



Slll-p— 



F (0,a,R) 




(16) 

Note that we pulled the a — ► asymptotics out in Eq. 
(fTS"]) . i.e. the functions Fi are finite for a = 0. For a and 
both nonzero, the integrals must be treated numerically. 

The form of the interaction Hamiltonian in (fT4|) (ob- 
tained by choosing a coordinate system where R = x) 
is convenient for reading off the various special cases 
of an RKKY interaction with no spin-orbit effects in- 
cluded (TJueis. ) ; the contribution from a pure Rashba in- 
teraction (i?R as hba)i the contribution from a pure Dres- 
selhaus interaction (Hj) Tess .), as well as the contribution 
to RKKY coming from the interference between the lat- 
ter two when these are simultaneously present (i?interf.)- 
While these expressions are suggestive for physical in- 
terpretations there is actually a more useful choice of 
coordinate system, obtained by choosing the angle be- 
tween R and y equal to — arctan(a//3). With this choice 
AR = a (0, (a 2 — (3 2 ) cosarctan(a//3)), and only t v ap- 
pears in the Green's function. The interaction then takes 
the simpler form 

#rkky = K u Si ■ S 2 + K lsing S v S% + Kdm (Si x S 2 ) y . 

(17) 

The three terms in (| 1 7|) can be identified with the Heisen- 
berg, Ising, and Dzyaloshinsky-Moriy a 37 ' 38 interactions, 
respectively. Their coefficients K^i, -Kising, -Kdm depend 
on the distance between the two spins, the Kondo cou- 
plings Ji,2, and the spin-orbit couplings a and 0. Since 
ffRKKY in (flT]) manifestly conserves U{\) spin symmetry 
(spin rotations around the y-axis), Eq. (TTTJ) is easier to 
work with than the expression in (|15p . where it is less 
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obvious how to exploit the U(l) symmetry. It is here no- 
table that the structure of the interaction in Eq. (jTTJ) is 
the same as that obtained for a pure Rashba spin-orbit 
interactions in Ref. 1291 . In particular, the interference be- 
tween the two types of spin-orbit interactions (ffmtcrf. ) in 
Eq. (fT5|) does not produce a new structure in (fT7|) . On 
the other hand, the parameters Kr, -ftTismg, -Kdm do de- 
pend differently on R, a and (3 in the general case. It 
is also worth pointing out that the special case |a| = \(3\ 
gives AR = 0, which means that only the Heisenberg 
term ~ Si • S2 appears in the RKKY interaction, like in 
the absence of any spin-orbit interaction. The restoration 
of the SU(2) symmetry for equal strengths of the Rashba 
and Dresselhaus couplings has been noted earlier—, and 
predicted to give rise to a persistent spin helix (PSH), a 
helical spin density wave of infinite lifetime. The emer- 
gence of a PSH was subsequently reported in an experi- 
ment on GaAs quantum wells^. 

In the following we will find it useful to express the 
RKKY interaction in terms of the parameters K v = 
±(ff H + ifi S i ng ) and e^K 1 - = Ki smg + iK DM as 

ffRKKY = RyStSl + \e ie K x {S z x + iSl){S z 2 - iSf) + h.c. 

Equivalently, we can express -ffRKKY in terms of a rotated 
second spin, 

S' 2 ^e^-S 2 e- ies K (19) 

in which case the RKKY interaction becomes^ 

ffRKKY = K^Si ■ S' 2 + (K y - K^)S\S' 2 V . (20) 

Note that in the special case of Eq. (fT6|) with only Rashba 
or Dresselhaus interaction present, the coefficient of the 
second term in (f20|) vanishes and the RKKY interac- 
tion in terms of a rotated second spin is purely of the 
Heisenberg-type. For our purposes, the RKKY interac- 
tion in the forms of (fT5|) and (|2H)) are most suitable and 
will be used in the following. 

III. KONDO EFFECT WITH SPIN-ORBIT 
INTERACTIONS 

Having analyzed the effect of spin-orbit interactions on 
the RKKY coupling between the two localized spins we 
we must now try to understand how spin-orbit interac- 
tions influence the Kondo interaction per se. The case of 
Rashba interaction was recently investigated in Ref. l3ll . 
and we here extend the analysis to include also the Dres- 
selhaus interaction. Since the non-trivial fixed point of 
the TIKM is unstable against various perturbations (in 
the language of the renormalization group) , we propose a 
slight modification of the setup that we have discussed so 
far. The modified setup is depicted in FIG. ([2]), and has 
the property that it protects the fixed points. In the 
absence of spin-orbit effects, and with the spinful dots 



operated in the Coulomb blockade regime (so as to sup- 
press charge transfer between dots and leads, as well as 
between the dots and the reservoir), we can write the 
effective low-energy Hamiltonian of the device as 

ff = ffkin + JSi ■ o-i + JS 2 ■ ct 2 + K(R)S 1 ■ S 2 . (21) 

Here er; = ip\ a (2^)7-^7/^(3^) is the electronic spin den- 
sity in lead i (i = 1,2, with t the vector of Pauli ma- 
trices). The Kondo coupling J is generated to second 
order in the tunneling rates, here assumed to be equal, 
J ~ Vi/Ud, with Ud the charging energy of the dots. 
Note that whereas in Eq. JT]) the spin-densities of the 
electrons in the central reservoir were coupled to both 
spins, here the spin densities in lead i couple only to 
the spin i. In effect, we have two single-impurity Kondo 
models coupled via the RKKY interaction K(R)Si ■ S 2 
between the impurity spins. The size of the RKKY cou- 
pling K(R) is assumed to be much larger than the Kondo 
temperature of the reservoir, hence the direct Kondo spin 
exchanges between the localized electrons on the dots and 
the electrons in the reservoir have been neglec ted in Eq. 
(|2Tjl . For more details on this, see Refs. flalla . 




Lead 1 Lead 



FIG. 2: (color online) The physical system that allows for a 
study of the non-trivial fixed point of the TIKM. The different 
V are tunneling rates. The dots are operated in the Coulomb 
blockade regime, where charge transfer between the reservoir 
and dots as well as between the leads and the dots is strongly 
suppressed. 

To find out about the interplay of Kondo, RKKY, and 
spin-orbit effects, let us begin by writing down the single- 
impurity Kondo model in two spatial dimensions (x, y) 
with an added spin-orbit interaction: 

ff = J d 2 k efcV>feVfe,<7 

+ j d 2 kj ' d 2 k' ^T'Vfc,,,,, ■ S (22) 

+ J d 2 kj d 2 k' (k,o-\H s _ \k',o-')^^ k >,«>, 

where (cf. Eqs. ©) and ©) 

ff s -o = (ak x - [3k y )T y + {f3k x - ak y ) T x . (23) 

Here \k,a) — ip^ |0) are momentum eigenstates with a 
the spin z-component, and t 1 are the Pauli matrices. As 
before, a and j3 are the couplings for the Rashba and 
Dresselhaus spin-orbit interactions, respectively, and we 
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assume a spherically symmetric free-electron dispersion 

The single-impurity Kondo model in the presence of 
only Rashba spin-orbit interaction (/3 = 0) was stud- 
ied in Ref. [3l| where it was found that the qualitative 
low-energy properties of the model are unchanged by the 
added interaction. This can be seen by writing the Kondo 
model with spin-orbit interactions in the form of a two- 
channel Kondo model, where the different channels of 
electrons have different couplings to the impurity spin. 
The spin-orbit interaction in this case is absorbed into 
the kinetic and Kondo terms of the Hamiltonian and no 
longer appear explicitly, ft is straightforward to apply 
the same procedure to the Dresselhaus interaction. 

As usual it is convenient to expand the electron fields 
in partial waves around the impurity site, which we here 
choose to be R = 0. Only the to = fields participate in 
the Kondo interaction (m € Z being the orbital quantum 
number), with the kinetic energy being diagonal in to. As 
shown in Appendix [C] the Dresselhaus- and Rashba- type 
interactions both couple each field to exactly one other 
field, i.e. a field labeled by m gets connected to a field 
with either m+lorm-1. This simplifies the picture 
considerably, as can be ascertained by the help of FIG. 

m 



X 

\ 

X 

^+2,1 



rf 



^-2,1 



^-2,T 
x 



Kondo 



Rashba 



Dresselhaus 



FIG. 3: (color online) The Kondo interaction as well as both 
kinds of spin-orbit interactions only couple electron fields with 
certain quantum numbers, as depicted here. In particular one 
sees that spin-orbit couplings preserve a "pseudospin" a with 
the values a =j, say, for the fields on the left-hand side and 
a =| for the ones on the right-hand side. Note that the Kondo 
interaction only couples fields of different pseudospin. 

Since there are no couplings in the kinetic or spin-orbit 
terms between fields on the left-hand side and right-hand 
side of the diagram in FIG. El it is possible to choose a 
basis of simultaneous eigenstates of the kinetic and spin- 
orbit Hamiltonian, where each state is a linear combi- 
nation of states from only the left- or right-hand side 
of the diagram, but not from both. For non-zero spin- 
orbit couplings, each of those fields contains one of the 
to = fields and thus couples to the impurity spin. If 



both kinds of spin-orbit interactions are present, there 
is an infinite number of fields coupling to the impurity. 
Unless there is some other mechanism that truncates the 
angular momentum states that may appear, it is not easy 
to make a statement about the resulting physics. In fact, 
an infinite number of fields coupling to the impurity is 
a similar situation as one would obtain for a long-range 
Kondo interaction in the absence of spin-orbit couplings. 

As long as only one of the spin-orbit interactions is 
present, the Hamiltonian of Eq. (|22)l can be rewritten as 



H 



(24) 



/,/' 



where / = ± is a flavor index, and a — L j is an SU(2) 
pseudospin- 1/2 of the electron fields, corresponding to 
the left- and right-hand sides of FIG. (O, respectively. 
The electrons are here expressed in terms of the energy 
E (for technical details see Appendix [Cj) . This model is 
known as the (anisotropic) two-channel Kondo models. 
In general, the Jf t f(E, E') interaction is neither diago- 
nal, nor has degenerate eigenvalues. A non-degenerate 
interaction (whether diagonal or not) is known to drive 
the two-channel Kondo model towards the single-channel 
Kondo model (plus one channel of free electrons) under 
renormalizationii.. The physical explanation for this is 
that if one channel of electrons couples more strongly to 
the impurity than the other, screening is fully achieved by 
that channel in the low-temperature limit. In this case, 
the other channel decouples and behaves like one of free 
electrons. Therefore, the effective low-energy model de- 
scribing single-channel Kondo exchange in the presence 
of spin-orbit interaction of either Dresselhaus or Rashba 
type, is the usual single-channel Kondo model (without 
any further interaction) , plus a channel of free electrons 
which decouples and may thus be dropped. 

The same argument can be applied to the TIKM. 
Physically, this is easiest to understand in terms of the 
nanoscale device introduced in the beginning of the sec- 
tion. Here, the TIKM is understood as two single- 
impurity models, coupled only via the RKKY interac- 
tion. In this picture it is clear that spin-orbit effects 
should not change the critical behavior. One can, for 
instance, start with two decoupled single-impurity mod- 
els, each of which flows to the single-channel fixed point. 
Then one adds the RKKY coupling, which does not af- 
fect how the leads couple to the impurities^. The TIKM 
in the presence of either type of spin-orbit interactions, 
which is equivalent to an anisotropic two-channel TIKM, 
can therefore still be described by the Hamiltonian in Eq. 
(|2I[) as far as the low-energy behavior is concerned. 
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IV. CRITIC ALITY IN THE PRESENCE OF 
SPIN-ORBIT INTERACTIONS 

Before discussing the effects of spin-orbit interactions 
at criticality of the TIKM, let us briefly review the most 
important results for the case when no spin-orbit inter- 
actions are present. In this case, an unstable fixed point 
has been found to exist for a particular value of the pa- 
rameters, i.e. A sa 2.2Tr-— . At this fixed point the 
impurity spins form a degenerate doublet between the 
antiferromagnetic singlet state |+— ) — | — h) and one of 

the triplet states |++) + | ) 44 , where we label all states 

in terms of eigenstates of the y-component of the spin 
operators, i.e. Sf\±)i = ±||±)j, i = 1,2. In general, the 
fixed point is unstable against breaking of particle-hole or 
parity symmetry 14 . If charge transfer between the dots 
is suppressed, the breaking of those symmetries become 
irrelevanlji*. Even with charge transfer present, breaking 
the spin SU(2) down to U(l) with an Ising-type interac- 
tion ~ S\S 2 produces only irrelevant operators, while a 
Dzyaloshinsky-Moriya type interaction allows a marginal 
operator 45 . The Dzyaloshinsky-Moriya interaction is not 
invariant under discrete rotations of ir around the x- or 
y-axis, thus this operator may appear in that case. 

With this, the stage is set to treat spin-orbit interac- 
tions around the fixed point. Since spin-orbit effects give 
rise to a two-channel model, which in turn renormalizcs 
to a single-channel model, and the RKKY interaction 
changes as derived in Sec. |TT1 the effective Hamiltonian 
can be written as: 

H =Hy_i u + J Si ■ CT\ + JS2 ■ &2 

+ K y S\S y 2 + ^e l8 K x (S^ + iSf)(£| - iS^) + h.c. 

^RKKY 

(25) 

As pointed out above, the Kondo screening behavior 
is not affected qualitatively; under renormalization 
the system flows to the same fixed point as without 
spin-orbit interactions. What may change is the value 
of the Kondo temperature Tk>2L. Since small deviations 
around the critical value of A = A 4 - = K y w 2.2T K 
are relevant, changing the Kondo temperature while 
keeping everything else fixed may drive the system 
away from the critical point 4 ^. To keep the system at 
the fixed point it may therefore be necessary to fine 
tune the interactions, which is in principle possible 
in the proposed nanoscale device. From here on, we 
assume any changes in Tk to be compensated by 
modifying A accordingly. The only way how spin-orbit 
interactions can then influence the critical behavior is 
by the symmetry-breaking SU(2) — > J7(l), by which 
new operators may appear. In terms of the RKKY 
interaction in Eq. JIgJ|, both (K y ^ K x ,6 = 0) and 
(0 ^ 0, A 4 - = K y ) break SU{2) down to 17(1). We will 
refer to the former as a longitudinal anisotropy, while 
we call the latter a transversal anisotropy. As we shall 



see, the effect of this symmetry breaking depends on the 
presence of charge transfer between the two channels 
of conduction electrons 4 ^. There are various cases to 
consider. We begin with the simplest case of purely 
transversal perturbations: 

a. Transversal anisotropics without charge transfer. 
If A 4 - = K y the RKKY interaction of Eq. <\%$ takes the 
simple form 

#RKKY = KSi ■ S' 2 . (26) 

We rotate the spins of the conduction electrons coupled 
to the second impurity spin to match the rotation of the 
impurity spin (Eq. (fT9| ) 

= (27) 

The kinetic energy is invariant under these transforma- 
tions, as are Si and It follows that the Hamiltonian 
takes the form 

H = H kiu + K{R)S 1 ■ S' 2 + JS 1 ■ a x + JS' 2 ■ ar' 2 . (28) 

Without the primes this is precisely the same Hamil- 
tonian as for the original (isotropic) TIKM (cf. Eq. 
(|21jl). This result does not mean that changing 9 is an 
irrelevant perturbation under which the system flows 
back to the isotropic fixed point. Instead, the fixed 
points for all values of 9 should be identified since they 
arise from the same Hamiltonian. It should be clear 
that a charge transfer term ~ is not invariant 

under this transformation, indicating that something 
different may happen in that case. We will return to 
this issue later, after considering the case of purely 
longitudinal anisotropics. Before proceeding, it may be 
worth noticing that only the singlet state is affected by 
the transformation above, turning it into |H — ) — -Mr\ — (-). 

It is easy to verify that this gives the expected spin-spin 
expectation value of (Si ■ S' 2 ) = — h. In the original basis 
this means that the unrotated spin-spin expectation 
value can take any value between — \ and +4. 

b. Longitudinal anisotropics with or without charge 
transfer. If K y ^ A 4 - and 9 = 0, the Hamiltonian 
cannot simply be reduced to the isotropic one, even in 
the absence of charge transfer between the two leads. 
However, this perturbation is known to be irrelevant 
from the boundary conformal field theory (BCFT) 
solution in Ref. [lj, even in the presence of charge 
transfer. In fact, breaking SU(2) to U(\) in this way 
does not produce a new leading irrelevant operator 
either, which means that for longitudinal anisotropics 
the system flows back to the isotropic fixed point under 
renormalization and the scaling behavior of thermody- 
namic quantities is unaffected. It should be noted that 
despite of this, changing K y while keeping A 4 - constant 
is a relevant perturbation. This is due to the fact that 
the shift (K ± ,K y ) -> (A 4 - + S, K y + 5), where S is 
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some small number, drives the system away from the 
fixed point K = K y w 2.2T^, and is hence a relevant 
perturbation 4 ^. In the linearized renormalization-group 
flow around the fixed point, this means that the irrel- 
evant direction must be perpendicular to this. Only 
K 1 -, K y -► K ± + S,K y - S can be irrelevant, i.e. K 1 - 
and K y need to be changed judiciously in order for the 
system to remain at the fixed point. 

c. Both longitudinal and transversal anisotropics, no 
charge transfer. As long as there is no charge transfer 
between the leads, the two previous results can be 
easily combined to determine what happens when both 
longitudinal and transversal anisotropies are present: 
By the unitary transformation, the Hamiltonian of such 
a system reduces to a Hamiltonian where there are only 
longitudinal anisotropies present. These are irrelevant 
as noted in the previous case. 

d. Both longitudinal and transversal anisotropies as 
well as charge transfer. As pointed out before, the pres- 
ence of charge transfer disallows the unitary transfor- 
mation that we used in the previous cases. Exploiting 
the BCFT solution in Ref. Q we find that an exactly 
marginal operator (Ji\hi<^ y in the BCFT language) is 
allowed when both transversal anisotropies and charge 
transfer are present. The isospin component h\hi of this 
operator is only allowed when the charges of the two 
channels are not separately conserved. The spin compo- 
nent <p y is allowed if the spin-symmetry is broken down 
to f/(l) and the symmetry under discrete rotations by 
an angle it around the x-axis and the z-axis is broken as 
well. This is only the case if 9 ^ 0, i.e. in the presence of 
transversal anisotropies. A charge transfer term between 
the two leads ~ ip\ip2 + V4V'1j along with the Kondo cou- 
pling between the impurity spins and the leads gives rise 
to an additional RKKY-interaction between the two im- 
purity spins, on top of the interaction mediated by the 
central reservoir. The marginal operator does not drive 
the system away from criticality, i.e. the RKKY cou- 
plings K v and K are unchanged (changes in their values 
are either relevant or irrelevant, see FIG [4]). Instead, we 
expect the phase 9 to be replaced by an effective phase 
6 c f[, made up of contributions by 9 and by the coupling 
to the charge transfer terms, which appears as J_ in[l41 
This is consistent with the fact, that the marginal oper- 
ator is allowed only for 9^0. The models for different 
values of 9 can therefore no longer be identified as be- 
fore. This suggests that there is now instead a line of 
fixed points parametrized by 9 e g, and connected via the 
marginal operator. 

e. Further away from the critical point. It should be 
kept in mind that the notion of an operator being irrele- 
vant only carries meaning in the neighborhood of a fixed 
point. For larger differences between K y and K the 
situation may be very different. In fact, for K 1 - = an 
entirely different quantum phase transition is expected 




FIG. 4: (color online) Qualitative RG flow of the anisotropic 
model. The solid dots and the solid line are known results for 
the isotropic model 4 ^. The gray area marks the system close 
to a different model of quantum dots coupled via an Ising in- 
teraction where different behavior is expected 4 ^. The dashed 
line which we argue for in the text separates the RKKY singlet 
from the Kondo screened regime. As an artifact of the scale 
at \K\ — > oo it appears curved to coincide with the screened 
fixed point (K ± — —K y ). Note that at both the singlet and 
the screened fixed point, the direction along the semicircle is 
irrelevant. We thus expect there to remain a finite separation 
(its scale being set by Tk) between the relevant flow towards 
the screened fixed point and the dashed flow towards the crit- 
ical point, as shown in the enlarged inset. The curvature of 
the dashed line is not meant to suggest any deeper knowledge 
about its properties; however, close to the isotropic (unstable) 
fixed point it follows the direction of irrelevant longitudinal 
anisotropies. 



to occur at a particular value of K y —. On the other 
hand, for K ^> and general K y < 0, the situation is 
very much the same as for K = K y : For sufficiently 
large values of K- 1 , the impurities form an RKKY sin- 
glet, while for K < —K y the system enters a phase 
where the two impurity spins are locked in an S = 1 state 
and get screened by two effective channels of conduction 
electrons 41 . This suggests that for any given K y < 0, 
there is a value K^. it where a phase transition similar 
to the K 1 - = K y case occurs. In fact, in the presence 
of particle-hole symmetry a phase transition must occur 
(see the discussion in Ref. LL4|, p. 9530). Furthermore, 
the same line of arguments as used in Ref. HH to solve 
the TIKM can be applied in this case (note however, that 
there it was supplemented by nRG and BCFT results). 
Since there is no other energy scale in the system, we 
expect that there is no additional phase transition away 
from K = 0. As illustrated in FIG.HJ this suggests that 
the irrelevant flow to the isotropic fixed point originates 
from an unstable fixed point at \K\ — > oo which separates 
the singlet RKKY and Kondo screened phases. 
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V. ENTANGLEMENT OF TWO 
RKKY-COUPLED SPINS IN THE PRESENCE OF 
SPIN-ORBIT INTERACTIONS 

Having explored how spin-orbit interactions influence 
the critical behavior of the TIKM, we now shift focus and 
turn to the question how these same spin-orbit interac- 
tions affect the entanglement between the two localized 
spins deep in the RKKY regime (where the competition 
from the direct Kondo interaction can be neglected) . As 
discussed in the introduction, this is an important prob- 
lem considering proposals 17 '— for using RKKY-coupled 
quantum dots as entangled qubits for quantum comput- 
ing. 

As a measure of entanglement we shall employ the 
standard concurrence of formation? 5 , from now on sim- 
ply referred to as concurrence. To obtain an expression 
for the concurrence in the present case it will be helpful to 
have access to the most general form of a reduced density 
matrix for a two-spin system compatible with conserva- 
tion of the y-component of the total spin. We can write 
it in terms of six real parameters 01, . . . , clq as 



a 1 \++){++\+a 2 \— 
a 3 \+-){+-\+a 4 \— 
(a 5 +ia 6 )\+-)(-+\ 
(a 5 -ia 6 ) |-+)(+-| 



)(- 
•)(- 



(29) 



Note that only five of these parameters are independent, 
due to normalization. We can therefore express the den- 
sity matrix in terms of the five expectation values 



(SfSI) - (S»S») 



3/ = 

D = (5f Sf) - (S?S%) 
3H = -(S 1 -S 2 ) 

x± = (S!)±(s$). 



(30) 



The way the parameters in Eq. (|30|) are defined allows for 
an easy reduction to various limiting cases. For example, 
H ^ 0,1 = D = X ± = corresponds to the SU(2) case 
with a pure Heisenberg interaction between the spins. 
Actually, in the system we consider, the parameters A ± 
are always zero by symmetry, i.e. there is no spontaneous 
magnetization, neither uniform nor staggered. Whereas 
in the derivation of our formula for the concurrence (App. 
IBl . we keep X ± in order to obtain a result that is valid 
also for more general situations, e.g. in the presence of a 
magnetic field, we shall drop A ± in the following so as 
to lighten the notation. 

We shall also find it convenient to work in the "magic 
basis" , given by 



|e 2 ) 
lea) 
|e 4 ) 



i|++>-i|— > 
i|+-)+i|-+) 



(31) 



In this basis the concurrence can be written as 



C = max {0, Ai — A2 — A3 — A4} . 



(32) 



where Ai are the decreasingly ordered square roots of the 
eigenvalues of p magic x (p mag i C )*, with p mag i c being the 
reduced density matrix in the magic basis^. After some 
algebra (see App. [B] for details) we end up with the 
following formula for the concurrence: 



C = 4 max {0,H + 21} 



(33) 



This implies that the concurrence can actually be ex- 
pressed in terms of a single parameter 



C = 



c{0,£}, 



where 



£ 



{{S\SD - (SlSl) - 3(SlS 2 y )) 



(34) 



(35) 



To calculate the concurrence of any state of the two- 
spin subsystem (arising from a total Hamiltonian con- 
serving U(l) spin symmetry) we can use the formula (|34|) 
derived above. Deep in the RKKY regime where the com- 
petition from the Kondo interaction can be neglected (the 
case discussed in Sec. UTJ) , the spins are described by a 
Hamiltonian of the form 

H = K y SlSl + ^e 16 'K ± (Sl+iSf)(SZ-iSZ) + h.c. (36) 

which has three possible ground states: a rotated singlet 
state 

tegletW = \ (|+-) - e i9 |-+» ((+-1 - e- l8 (-+\) , 

(37) 

an Ising ground state 



and a mixture of rotated triplet states 



(38) 



PtripletW = Ringlet ( d ) + \ (++I + \ I > ( 1 . (39) 

corresponding to the cases K > —K v , < —K y (or 
\K y \ > K ± =0) and K 1 - = -K v , respectively. This can 
be read off immediately from the level scheme in FIG. [5] 
Note that the states in Eqs. (f3"T)) and (f3"9"]) are just the 
conventional singlet (9 = 0) and triplet states (9 = tt) 
in a basis where one of the spins has been rotated with 
respect to the other (cf. Eq. ([2H|) '). The parameter £ 
that determines the concurrence in those three cases is 
(cf. Eq. 



"^singlet 
^Ising 
rip let 



-l 

-1/3 
1. 



(40) 
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l+-) + e ifl |-+> 



of Sec. |TV] we find (cf. Eq. (J35J) 



!+-),!-+) 



!+-).!-+). I— ). 



I++>.I~ ) 



-K y 



K 



\++),\—) 



FIG. 5: Level splitting due to the RKKY interaction. 



hence the concurrence is unity for the singlet while it 
always vanishes in the latter two cases (see Appendix [B] 
for details). 

To understand what effects a spin-orbit interaction my 
have on the entanglement in the RKKY regime, consider 
the situation described by the iVamphtude in Eq. (Til)]) , 
i.e. an RKKY interaction given by 



H = -^R2 ! W S ™ 2R \l k F 



h 4 



-Si ■ S 2 . (41) 



It is clear from this that a change in a may change 
the overall sign of the interaction and thus drive a sys- 
tem with a maximally entangled singlet ground state to 
one with a non-entangled triplet ground state, or vice 
versa. As an illustration, one may consider a double- 
quantum dot structure patterned in an InAs heterostruc- 
ture, which, due to its large electron mean-free path, is 
favored in most spintronics applications 4 ^. Using data 
for a heterostructure grown by molecular beam epitaxy^, 
with a gate controlled Rashba parameter tunable in the 
range 1 — 5 x 10~ n eV m, and with an effective mass 
?7i ps 0.4m e , one finds that the inverse Rashba spin-orbit 
length ma/ti 2 varies in the interval 0.01 — 0.1 A -1 (to 
be compared to the Fermi wave number kp ~ 0.02 A -1 , 
and the distance R between the dots, say, of the order of 
magnitude R ~ 10 2 A). In the more general case — when 
the F\- and ^-amplitudes in Eq. (fT6| come into play — 
it must be expected that a transition to a system with a 
non-entangled Ising ground state can also occur. We con- 
clude that spin-orbit interactions may drastically reduce, 
or, enhance, the entanglement between two spin qubits, 
with the precise effect depending on materials and design 
parameters, as well as on the strength of applied electrical 
fields. It is here important to emphasize that estimating 
the value of the RKKY coupling in a real device is no 
easy task. In addition to the uncertainty in determining 
the distance between the two spin qubits, there are diffi- 
culties in judging the influence from the mean-free path, 
dephasing length, etc. An exquisitely precise control of 
the experimental set-up is a sine qua non for a reliable 
estimate of the RKKY coupling on which entanglement 
depends so crucially. 



VI. ENTANGLEMENT AT CRITIC ALITY 

Using the formalism developed in the previous section 
(see Appendix iBl for details), and exploiting the results 



e, 



critical 



(42) 



i.e. the concurrence vanishes at the critical point K = 
K y pa 2.2Tk as in the 5/7(2) symmetric case with no 
spin-orbit interaction 2 ^. It should be noted that in the 
cases of a vanishing concurrence in the noncritical RKKY 
regime, i.e. for the Ising doublet and the triplet states 
(discussed in Sec. IV), the value of £ was always a fi- 
nite and negative number, and the concurrence vanished 
by virtue of taking the maximum of zero and this num- 
ber. In the case of the critical ground state, £ is by itself 
zero, which means that a small deviation from the criti- 
cal point could lead to a non-zero value for the concur- 
rence. In other words, the point where the concurrence 
first vanishes marks the critical point. Here a stronger 
antifcrromagnetic coupling drives the system to the free 
fixed point where the spins form a fully entangled sin- 
glet, while a weaker antiferromagnetic or even ferromag- 
netic coupling drives the system to the fully screened 
two-channel Kondo fixed point where there is no entan- 
glement between the spins. 

It is important to keep in mind that the BCFT result in 
Ref. fil about the irrelevance of anisotropies is only valid 
for small deviations around the isotropic case. In the de- 
vice which we suggest (cf. FIG. 3), we have shown that 
the anisotropy in the sz-plane can be arbitrarily large, 
i.e. the phase of K does not change the picture. For 
large anisotropies along the y-axis this will be quite dif- 
ferent. We are not able to describe the exact behavior in 
the entire parameter plane spanned by K and K v , but 
for certain regions there are strong arguments that the 
system should behave in a certain way. Let us first con- 
sider the case \K V \ -c K . For very large values of K 
the impurity spins will form a singlet and the conduction 
electrons are decoupled, just as for K v ps K ^> Tjc- 
For K -> (and \K V \ < K^) the system will flow to 
the two-channel Kondo fixed point, again just like in the 
isotropic case. In fact, the same line of arguments used 
in Ref. |H can be applied. The impurity state at the crit- 
ical point should be the same, in particular also having 
vanishing concurrence. We expect that this fixed point 
can be connected to the isotropic fixed point by tuning 
both K y and K appropriately, and extended for larger 
values of K 1 - > -K« > 0. For \K y \ > K ± , \K y \ > T K 
we are in the RKKY regime where the impurity spins 
are completely decoupled from the conduction electrons. 
We thus know that K y <C gives vanishing concurrence, 
while for K y >0a nonzero value of K 1 - raises the con- 
currence from to 1. There is, however, an intermediate 
region K v ~ TkjK 1 - = where a different quantum 
phase transition is expected to take placed. An analy- 
sis of the crossover to this quantum critical point, or of 
the entanglement properties in its vicinity, is beyond the 
scope of this paper. 
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VII. SUMMARY 



The RKKY interaction is proportional to 



We have carried out a systematic study of how spin- 
orbit interactions — of Rashba as well as of Dresselhaus 
type — influence the Kondo and RKKY interactions in 
a nanoscale device described by a two-dimensional two- 
impurity Kondo model. By absorbing the Rashba and 
Dresselhaus interactions as anisotropies of an effective 
RKKY coupling we can monitor their effects for varying 
coupling strengths. This provides us with a controlled in- 
route to study the quantum critical behavior of the model 
in presence of spin-orbit interactions, exploiting known 
results from boundary conformal field theory^ and an 
effective bosonized Hamiltonian approach 44 . Most strik- 
ingly, and as illustrated in Fig. 4, for a particular com- 
bination of symmetry breaking terms there is strong evi- 
dence for a line of critical points exhibiting the same uni- 
versal behavior as that of the known isotropic model (up 
to RG irrelevant scaling corrections). It remains a chal- 
lenge to rigorously establish the existence of this critical 
line. We have also studied how spin-orbit interactions 
influence the entanglement between two RKKY-coupled 
spinful quantum dots (alias "qubits" in the language of 
quantum computing). Using data for a nanoscale device 
patterned in a gated InAs heterostructure we find that a 
gate-controlled Rashba spin-orbit interaction may drive 
a maximally entangled state to one with vanishing en- 
tanglement, or vice versa (as measured by the two-qubit 
concurrence). This has important implications for pro- 
posals using RKKY interactions for nonlocal control of 
qubit entanglement in semiconductor heterostructures. 
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Tr [(Si • t)G(R)(S 2 ■ t)G(-R)} . (A3) 

Inserting G from Eq. (|A2|) into (|A3|) and using the trace 
formulas given above, we obtain: 

Tr[(SyT)G(R)(S 2 ■ t)G(-R)] = 

2GqSi ■ S 2 
- AGoGna (Si x S 2 ) y 
+ 4G GnP {Si x S 2 ) x (A4) 
+ 2G\a 2 (Sx ■ S 2 - 2SfSf) 
+ 2G\f3 2 (Si ■ S 2 - 2Sf Sf ) 
+ AGjaP (Sf S v 2 + SfS% ) . 



APPENDIX B: CONCURRENCE FOR VARIOUS 
STATES 

The concurrence for two qubits (two-spin subsystem) 
can be written as 



C = max {0, Ai — A2 — A3 — A4} . 



(Bl) 



where Ai are the decreasingly ordered square roots of 
the eigenvalues of p ma g ic x p* nagic , with p magic being the 
reduced density matrix in the magic basis^. In terms of 
H, I, D, and X defined in Eq. (|3"01) . it is given by 



'-H' - 21 -1X+ 

iX+ —H' - 21 

4/ 

D - 






- H' D + iX- 
-iX- \ + 3H 

(B2) 

where H' = H — 1/4. Since /3 mag ic X Pmagic i s block 
diagonal it is straightforward to calculate the eigenvalues 
and one finds 



APPENDIX A: RKKY INTERACTION FOR R = x 

To derive the form of the RKKY interaction, it is neces- 
sary to perform traces over various combinations of Pauli 
matrices t x ,t v , and t z . Specifically, we need the rela- 
tions 

Tr [(t ■ A) (r • B)] = 2 A ■ B, 
Tr \(t ■ A)(t ■ B) r'j = 2i(A x B) u 

Tr \{t ■ A) r j (t ■ B) t'J = 2A z Bj +2A j B i ~28i J A- B, 

which can be easily verified using 

(r • A) (t ■ B) = (A- B) +i(A x B) ■ r. (Al) 

For the case R = x, discussed in the text after Eq. (fTB")) . 
the single-electron Green's function is 



\2 
A l,2 



H -21 



(^ 



D 2 - (X- 



(2I-2H) 



21 + H 



±2 



\ 



((2/ — 2H) 2 + D 2 





Note that the parameters are not fully independent; all 
eigenvalues are positive numbers. 

To calculate the concurrence by the formula in (|F31|) . 
it is necessary to order the eigenvalues by size. Since A_ 
can never be the largest eigenvalue, there are only two 
cases, A+ > Ai and Ai > A+, to analyze. First consider 
A+ > Ai. In this case Eq. (|32[) takes the form 



G = G + G ia T y -Gi!3T x 



(A2) 



C = max{0, A+ - A_ - 2Ai}, 



(B4) 
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which, in particular, implies that the concurrence van- 
ishes for A + — A_ < 2Ai. As for the second possibility, 
Ai > A + , inspection of Eq. (|B1|) reveals that the concur- 
rence vanishes identically for this case since Ai — A2 = 
and A± > 0. However, A + — A_ < A + < Ai < 2Ai, thus 
the concurrence also vanishes by Eq. (IB4|) . It follows 
that the second case is already contained in that formula 
and it is not necessary to treat the two cases separately. 

In the following we shall specialize to the case with 
X ± = (cf. Section V) in which case the Ai simplify to 



Ai = - - H - 21 



A± = ^D 2 + (2I-2Hy± 



21 + H 



(B5) 



and we immediately end up with Eq. (|34[) , i-c 



C = max < 



47- 



2H + \ 



-1+2H- 



41 



4max{0, J ff + 21} 
max {0, £ } . 



(B6) 



To calculate the concurrence for a given reduced den- 
sity matrix in the magic basis, as parameterized in Eq. 
(IB2|) , we need to calculate the expectation value £ as de- 
fined in Eq. 1351 The corresponding operator takes the 
following form in the magic basis: 



diag(-l,-5,3,3) 



(B7) 



Deep in the RKKY regime there are only three possi- 
ble states for the two coupled qubits:, a singlet, an Ising 
doublet, and a triplet. Starting with the singlet state, 
(Eq. [37]) . its density matrix in the magic basis is given 
by 



"let 



(9) = diag 








sin z § i sin 9 
isin(9 cos 2 I 



(B8) 



and £ is determined as 



"^singlet (0) = Ti[O s 

Psinglct 



1 

-Tr 

3 



2 



3 sin' 
I sin ( 



§ sin0 
3 cos 2 f 



= 1. 



(B9) 



This means that the concurrence for the singlet state is 
unity, independent of 9. In other words, the well-known 
fact that a singlet state is maximally entangled is imme- 
diately reproduced in our formalism. The Ising doublet, 
(Eq. |3"8|) , has the density matrix 



By inspection £ 



Plsing " 
1 



diag(l/2,l/2,0,0). 



(BIO) 



T, thus the concurrence is 



parameters for ptripiot(^) can easily be obtained from the 
ones for p s ingiet(#) and pising, since these density matrices 
as well as the spin-operators are block diagonal. One 
finds 



£triplct(#) — 



"-let 



(0) 



If 

g rising — g , 



and the concurrence is again seen to vanish identically. 

Turning to the critical ground state of the TIKM, its 
reduced density matrix for the two qubits is given by 
Eq. (|B11[) . with a non-zero value of 9 when spin-orbit 
interactions are present, 



Pcritical(^) — J 



(I 
1 






sin( 






2 sin 2 ^ 
\0 sin 9" 2 cos 2 !/ 



(Bll) 



From this we find that ^critical = 0, i.e. the concurrence 
vanishes at the critical point. 



APPENDIX C: 2D KONDO MODEL IN THE 
PRESENCE OF RASHBA OR DRESSELHAUS 
INTERACTION 

The Hamiltonian for Rashba and Dresselhaus spin- 
orbit interactions in two dimensions, which we take as 
the x-y plane, is (see Eqs. © and ([5])) 



H S . Q = (J3k x - ak y ) t x - (fiky - ak x ) ■ 



(CI) 



We can define raising and lowering operators for spin 
and angular momentum as 



r ± = ~(t x -it v ) 
2 V ' 

L = k x it iky , 



to rewrite the Hamiltonian as 

H s . = ar-L+ + f3r+ L+ + h.c. 



(C2) 



(C3) 



This expression suggests to write the free fermion fields 
tph in terms of an angular momentum quantum number 
m and the magnitude k of the momentum as 31 



2 71 



(C4) 



The second-quantized spin-orbit Hamiltonian can then 
be expressed as 

H s . = [dk T k 2 (a^ m+1 (r + rJ^ m ^ 



^ViCOrVwO+kc. (C5) 



zero for this case. As for the triplet state, Eq. (13^1) . the It follows that the 2D Hamiltonian for the Kondo model 
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with added Rashba and Dresselhaus spin-orbit interac- 
tions takes the form 



H = H kin + H K , 



ondo 



Hs-o 



(C6) 



where 



#Kondo = J I dk 



dk'Vkk 1 ^2 ipl a o T « Vv,o,< 



s. 

(C7) 

and where -ff s _ is given in Eq. (|C5|) . 

As manifest in Eq. (|C5J) . when both kinds of spin-orbit 
interactions are present, all m-states couple to the impu- 
rity. On the other hand, if only one type of spin-orbit 
interaction is present, there are just two states coupling 
to the impurity. To diagonalize the spin-orbit Hamilto- 
nian in this case, we define new fields as 



1pk,m,l,± = (V*a,m ± ^ {n ~ 1),2 i>kA,m-7^ , 



(C8) 



where n = 1 in the Rashba case ((3 — 0) and n = — 1 in 
the Dresselhaus case (a = 0), and where the f,J, labels 
on the new fields refer to the "pseudospin" a introduced 
in FIG. (4). 

In terms of these fields, the three terms of the Hamil- 
tonian in (|C6p are given by 



-Hkin = 



#Kondo 



k,m,f' l r > k,m,(T,fi 



/,/'=± 



m ct=TJ 



(C9) 



where = a if the spin-orbit interaction is of the Rashba 
type and /x = j3 if it is of the Dresselhaus type. Since all 
terms are diagonal in m and only the m = fields take 
part in the Kondo interaction, we drop all other fields 
and suppress the m quantum number which is always 
zero from here on. Defining = e& ± fj,k, we combine 
-ffkin and H s . to 



#kin = / dk ^ 4^kJ^k,S,f- 



(CIO) 



As a last step we define a set of fields as 
J dkJl6(E - e{)^ 9j 



EJ,S 



J dk'6(E 



(Cll) 



to obtain the Hamiltonian 



H = 



E,f,u 



/,/'=± 



~ fdE f dE' J fJ '(E,E')^i f T^ E ,j,^-S, 



/,/'=± 



(C12) 



where 



J fJ '(E,E') = J J dk6{E-e f k ) J dk'5(E'- 



-<£')■ ( C13 ) 
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